The shape and mechanics of curved fold origami structures 
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Abstract - We develop recursion equations to describe the three-dimensional shape of a sheet 
upon which a series of concentric curved folds have been inscribed. In the case of no stretching 
outside the fold, the three-dimensional shape of a single fold prescribes the shape of the entire 
origami structure. To better explore these structures, we derive continuum equations, valid in 
the limit of vanishing spacing between folds, to describe the smooth surface intersecting all the 
mountain folds. We find that this surface has negative Gaussian curvature with magnitude equal 
to the square of the fold's torsion. A series of open folds with constant fold angle generate a 
helicoid. 



Introduction. — In the hands of an artist, a sin- 
gle sheet of paper can be folded into a seemingly end- 
less variety of three-dimensional, and strikingly beautiful, 
forms. Among the various fold patterns, structures built 
from curved folds, first folded in the late 1920's by a stu- 
dent at the Bauhaus [JL^, stand out as being particularly 
stunning [2]-[5]. The mechanics of these structures arises, 
at least in part, from the interaction between the avoid- 
ance of in-plane stretching within the sheet and the three- 
dimensional geometry [6]. Besides its own artistic value, 
origami has, for many years, shown itself to be a valuable 
tool to approach problems in architecture, structural me- 
chanics, and engineering [7HTT]. This new framework of 
curved pleated structures has already been demonstrated 
to be very useful for designing materials with new mechan- 
ical properties [TT] (e.g. anisotropic response to external 
forces and structural stiffness). These techniques could 
also potentially serve as a platform for design at many 
different scales. 

In this article, we will explore the shapes adopted by 
multiple, circular curved folds such as those shown in Fig. 
[H In our analysis, we shall assume that the sheets from 
which the folds are made resist stretching between the 
creases, which implies that their shapes can be approxi- 
mated to developable surfaces (surfaces with zero Gaus- 
sian curvature |l2]). Since it is not clear that this is always 




Fig. 1: A curved- fold origami structure with three prescribed 
concentric circles that is either (a) closed or (b) open (non- 
constant curvatures). The top figures demonstrate computa- 
tional origami, based on three concentric, circular folds of fixed 
spacing, wHi]^ = 0.1. (a): the inner fold has curvature 2K^gl\f?) 
and torsion t^(^s) = 1.26sin(2s), chosen so that the origami is 
closed, (b): the initial curvature is k^{s) = 1.2^c^ and torsion 
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possible (a fold made from nested squares is not foldable in 
this strict mathematical sense [l3]), we explore the possi- 
ble presence of singularities within the sheets. In fact, our 
findings suggest that a sufficiently narrow spacing between 



p-i 



M. A. Dias et al 



folds is generically sufficient to prevent singularities from 
appearing. We show that nested helices are a solution for 
open concentric circular folds, such as in Fig. [TJd. Fi- 
nally we consider the asymptotic shape of the surface that 
passes through all the mountain (or valley) folds in the 
limit that the number of folds diverges. The nested helix 
solution, in this limit, becomes a portion of a helicoid. 

Geometry of multiple folds. — To establish our no- 
tation, we review the geometry of a single curved fold [HI], 
and extend these results to recursively construct multiple 
pleated structures. Let Ci{s) be the embedding of a set of 
folds, labeled by i, in space. We choose s to be the arc 
length of the first fold, and define 5* to be the arc length 
of the i^^ fold. Consider the diagram in Fig. [2^ showing 
three consecutive folds. The tangent planes on either side 
(dz) of the i^^ fold, where (+) and (— ) respectively denote 
the outside and the inside sheets (Fig. [2)3), are spanned 
by the unit tangent V = dgid" and its orthogonal unit vec- 
tors Uj^. For completeness, we also define N!^ = t* x u!^ 
to be normals to each tangent plane. Since we assume 
the surfaces on either side of the i^^ fold are developable, 
they are spanned by two unit-length generators respec- 
tively. The relationship between two consecutive folds is, 
therefore. 



c\s) = e-\s)+vX\s)^^\s), 



(1) 



where 



gi {s)=i' [s) COS [7i [s)] T u'± [s) sin [7^ (s)] , (2) 

where the variables 7^ are the angles between the genera- 
tors and unit tangent. The distance v^_^^{s) is the distance 

^ fold 



along the generator from the {i — 1)^ fold to the 
so it depends implicitly on 7!^. 

It can be seen from the diagram in Fig. [2)3 that the 
decomposition of the vectors u!j-(s) are as follows: 



Vi±{s) = n*(s) sin 



± h\s) cos 



, (3) 



where n* = dgiV /\dsit'^\ and b* = t* x n* are the Frenet- 
Serret normal and binormal, respectively, and 0^{s) is the 
dihedral angle of the fold. The triad of vectors {t% n*, b*}. 
Fig. [2)3, forms natural moving frames on each i^^ fold, 
where curvature and torsion are respectively defined as 
follows, n\s) = h\^s^V and t'{s) = -h\^s^h\ The 
angles 7^ can be expressed in terms of geodesic torsion. 



Tg^{s) = r* ^ and normal curvature, 

^tv^ COS [0^2] , as 



cot = ±^ 



is) 



(4) 




Fig. 2: (a): A schematic representation of the geometry of 
the generators in the multiple folded configuration, (b): A 
mountain fold split along the crease in order to demonstrate 
the Frenet-Serret frame, which is the same for either side, and 
the Darboux frames for the regions (+), on the left, and (— ), 
on the right. The generators are also shown. 



77^ (s) = cot [7^(5)], such that 



4 cot [9^{s)/2] 



(5) 



relates any two sides, (+) and (— ), across a fold. Yet 
another curvature that needs a definition is the geodesic 
curvature, n^g = tz^ sin [0'^/2]. This quantity is the pro- 
jection of the space curvature, onto the tangent plane 
and can be interpreted as the curvature of the crease pat- 
tern before folding [4J . It is important to notice that is 
constant in s\ for every i, if the crease pattern is made of 
concentric circles, similar to the patterns that generated 
the structures in Fig. [H 

Fig. [2^ shows the schematics of the fold geometry, 
which is constructed so that the generators of consecutive 
curves are related to each other by 



(6) 



When the crease pattern is given by a sequence of concen- 
tric circles evenly spaced with separation simple planar 
geometry allows us to write an expression for v'^^ for the 
distance to the {i ± 1)^^ fold from the z^^. 



= ^sin7^ ± ±2) +sin^7^, (7) 

where we have suppressed the s dependence for brevity. 
Finally, inextensibility requires 



(8) 
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From these definitions, it follows that r*^ = r*_ — dgiO'^ 
It is convenient to define the variable 



and K,' 



7V+ 



N- 



which resembles our choice to prescribe folds such that 
These equations provide a means 
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of starting from a single fold, say at c^, and generating 
the shape of all subsequent folds. First, we define the 
parameter li{s) such that dgi = ds/l'^. Then requiring 
that dgic'^ be normalized yields 



= (1 




>rV\/i+(^r)' 



i+(r?r 



(9) 



Using the recursion relations for the embedding of the 
curves, equation ([T]), and the explicitly form of the rate 
of rotation of the frame {t%u!^,N!^} (Fig. [2}3), 
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we can derive the following relationships 



and 



t'~^ 



(11) 



(12) 



Then, the equation (j4]) yields 



Finally, we have that 6^ = 2 cot ^ (^i^N- I 

Eqs. dHJ, (pT]) . and (p!2]) can be solved numerically for 
either closed or open fold patterns, as shown in Fig. [T^-b. 
However, each subsequent fold requires higher numerical 
precision, so that structures with more than five folds be- 
come prohibitively difficult to evaluate computationally. 

It has been an open question whether closed folds of 
the type shown in Fig. [1^1 can exist without stretching 
beyond a single fold [13 . The numerical evidence depicted 
in Fig. [1^1 suggests the existence of structures, at least for 
a finite number of consecutive folds. To go beyond this, 
one must understand what kinds of obstructions can occur. 
Since the surfaces between folds are developable [12] , these 
obstructions arise as singularities because the generators 
cross. At these singularities, an elastic sheet will stretch at 
great energetic cost to avoid a diverging bending energy. 
In our notation, they can be expressed as bounds on how 
quickly the 7^ can change [6]: 



(14a) 
(14b) 



Additional constraints arise because generators on the in- 
side sheet of each fold may emerge from the crease at a 
value of 7L sufficiently small that the generator fails to 
meet the previous fold. This geometric constraint, which 
occurs when the discriminant in equation ([7]) is negative, 
can be translated into a bound for the torsion. 
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: cot 



2wK,i 



(15) 



Mechanics of Helical folds. — We can find solu- 
tions to the recursion relations for open folds of constant 
curvature and torsion. Because constant curvature im- 
plies a constant dihedral angle, the equation ^ becomes 
T]!^ = rf_' The recursion relations (j9]), (pTj) . and ([12]) then 
simplify to the following relations 



(i + ^4)r-^ 



(1 
1 



Tn = 



1 



■ WK 



(16a) 
(16b) 
(16c) 



In Fig. [TJd we show a solution for three consecutive folds 
generated by an initial open helix. A singularity occurs 
on the n^^ fold when 



1 — w/{r -\- nw) 



(13) \/{r^nw)\Tn\ = -==== 



cot 



(17) 



As long as these bounds are always satisfied, there is a so- 
lution to the recursion relations corresponding to a curved 
fold. In fact, it appears that the constraints can always 
be satisfied for sufficiently small w. Since vl^'^ w/ sin 7^ 
ioT n^gW <^ 1^ w can be decreased so that the right hand 
sides of Eqs. (p!^ and ([15]) become arbitrarily large on 
folds adjacent to the i^^ fold. Whether the structure must 
terminate after a finite number of folds for any particular 
spacing of course, is more difficult to determine. 

We formulate the mechanics of these structures by 
adding the individual contributions of bending energies 
coming from each pleat as well as fold energies coming 
from each crease i. The developable sheets on either side 
of each fold have bending energy [6] 



'b± 



sin 7^ 



(18) 



X In 



sin 7^ - vi^ {f^i T7±) 



where B is the bending stiffness. For a multiply folded 
structure, we perform a constrained sum over adjacent 
pairs of folds in order to avoid double counting facets when 
summing over i. Therefore, the energy for a structure with 
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Fig. 3: Energy landscapes normalized by B. We have set — 1 and Oq — 27r/3 for alH G {1, 2, n}. The black dots on each 
energy landscape represents the energy minima at t^/a^J =0 and / n^g = 1.15, 1.14, 1.13, 1.12, 1.11, and 1.1 respectively for 
= 0.1, 0.2, 0.3, 0.4, 0.5, and 0.6. The red dashed lines mark the limit of the singular region determined by Eqn. (|15|) . 



n folds is given by 

( (n+l)/2 

E 

i=i ± 

(n-l)/2 



Eb = ( 



As in [6l, we use 
2 



Er. 



E E^, 

=1 ± 



dsV 



for n odd 



^ for n even. 



(19) 



cosl- 



cosl- 



(20) 



for the fold energies, where <^ is a measure of the stiffness 
of the fold and 9q is the preferred angle for each crease. 

Typical energy landscapes are shown in Fig. [3l These 
were calculated for three fold structures, like the one 
shown in Fig. [TJd using the formulas (fT9|) and ([20]) . Sin- 
gularities appear along the dashed curve in Fig. O as 
given by the relation ([15]). Beyond this limit, stretching 
becomes necessary in order to avoid a diverging bending 
energy. Finally, the minima of energy all have =0, as 
can be verified by folding paper. 

Continuum Limit. — Artistically, origami is used 
as a medium for sculpture. In that sense, one feels a 
distinct impression that a folded sheet approximates a 
smooth shape. An origami structure made from many 
nested squares, for example, has the appearance of a sad- 
dle and, in fact, asymptotically approaches a hyperbolic 
paraboloid (when extra creases are added) as the number 
of folds becomes large fT3^. A structure with concentric, 
circular folds also has the impression of a discrete approx- 
imation of a smooth saddle (Fig. (1^) and, in this section, 
we seek equations to describe the resulting smooth shape. 

To be more precise, consider a smooth surface passing 
through all of the mountain folds of an origami sculpture, 
as shown in figure Fig. |4l We seek equations governing the 
shape of this surface as the width between folds decreases 



but the number of folds diverges. To proceed, we consider 
two mountain folds and the valley fold lying between them. 
Since w is small, we only express results up to O {w). Then 



2w 



2/^^^cot (i9V2) 



t*-sin|^-Jn* + 0(^). 

(21) 

Taking the limit where w ^ 0^ we let i become a contin- 
uous independent variable, denoted x- This yields a con- 
tinuum approximation for the origami structure that can 
be interpreted as the initial mountain (or valley) crease 
evolving with a "velocity" defined by 



V(5, x) = d^cis, x) = (^{s, x)t(s, x) + X)n(s, x), 

(22) 

where the collection of all "snapshots" along the variable 
X traces out a continuous surface. Here we define the 
components of V(s, x) as 



<^(5,x) 
P{s,x) 



dsO{s,x) 



2l{s,x)^9{x) cot [0{s,x)m 
~0{s,xy 

sm 



(23) 
(24) 



These components can be entirely expressed in terms of 
the curvatures through the constraint = sin~^ {tig/ n). 
Since the discrete variation ([2T]) and the velocity ([22]) are 
spanned by {t,n}, we note that they both lie on the os- 
culating plane of each curve. Therefore, the osculating 
plane becomes the tangent plane to the approximate sur- 
face when the limit ^ is taken. 

The above discrete equations, ([8]), ©, ([U]), and ([T2]) . 
turn into continuous evolution equations for i^g{x)i K^^X)i 
n{s^x) and r(s, x) of the mountain folds in the continuum 
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hmit, 



--dsT ^2--dsP 



d^Kg = —n^g (25a) 

d^l = dsa - Inp (25b) 

-T^)p^^dsn (25c) 

2KT(3-jdsr. (25d) 



The above set of equations (|25]) . with appropriate initial 
data, forms a complete set of PDE's. Before looking for 
specific solutions, we establish some general formulas for 
the geometrical quantities on the smooth, approximating 
surface. We compute the surface's first fundamental form 
to be 



I = ds'^ + 2la dsdx + (a^ + /3^) ^ix^ 
and second fundamental form 



II = -21 /3t dsdx ■ 



Thus, the mean curvature is 



Ik 



21k? 



and Gaussian curvature 



K = 



(26) 

(27) 
(28) 

(29) 



This final result implies that < 0, with K = {) only for 
a fiat surface. 

When n and r are independent of s, d^l = l/x so that 
^ = X-> where we have chosen s to be arc length when 
X = 1. Thus, 



X 



so that T = To/x^- Finally, this implies that 



2 2 
— T 



(30) 



(31) 




Fig. 4: Helicoid surface approximating helical pleated struc- 
ture. The translucent surfaces are the solutions of the con- 
tinuum limit c{s,w), for kq = 2 and to = 1. The top (blue) 
surface touches the mountains of the pleated surface, while the 
bottom surface (red) touches the valleys. 



The helicoid solution suggests that, in order to obtain 
the asymptotic behavior of more complex structures, one 
should use the following rescalings for /(s,x), f^i^^x) 
r(5,x)* I = Vx? ^ = X^ f = rx^. Thus, we have 
a = jds{0/2) tdin{0/2) and the evolution equations 



d^K = 



= -dsa 
X 

fdsR 



I 



_J_dsf 

X^^ I x^^^^ 
adsi^ 1 r ^ 
X I X ^ 



(33a) 
(33b) 
(33c) 



As X becomes large, these equations reduce to the asymp- 
totic equations 



This equation has solution hz = Vv^o+^o)x^^o /x^ • 
Notice that H = and K = — = — Tq/x^ varies only 
along X- This information is already sufficient to deter- 
mine that the shape must be a member of the helicoid- 
catenoid family. In fact, it can be verified that the surface 
embedding is given by 



c(s,x) 



sin (^s^/K^~Tr^ 



(32) 



which is indeed a helicoid. Fig. |4] shows how well the 
solutions ([32]) approximate the corrugated surface. 



l = X 

T =To{s)/x^ 
,^-dl{l/-K). 



(34a) 
(34b) 
(34c) 



Writing F = 1/^, we see that this last equation is a non- 
linear diffusion equation of the form d^T = T^d^V. In the 
large x hmit, F will become x- independent, implying that 
there is a n^i^s) such that it satisfies —dsi^o/i^o = C asymp- 
totically, where C is a constant number. The only periodic 
solution is a^o equal to a constant. This means that ini- 
tial curvatures fiow towards a fixed point, the constant 
value tZQ. Interestingly, this scaling form for the curvature 
implies that becomes nearly constant in s and x as x 
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becomes large. This satisfying result indicates that con- 
centric, circular folds are not overly frustrated from the 
point of view of their fold angles. In particular, we expect 
9 to be very close to its prescribed value, even when the 
fold angle cannot be exactly constant because of geomet- 
rical constraints. The mean curvature of the asymptotic 
surface is 

1 dsTQ^s) 



H{s) = 



(35) 



Conclusion. — In this article we explored the ques- 
tion of designing pleated folds. We were particularly in- 
terested in concentric circles crease patterns which are 
folded by alternating mountains and valleys. We pro- 
vided the construction of two examples, closed folds that 
yield a saddle-like shape and open ones resulting in helical 
shapes. Closed folds are particularly interesting objects 
because once they are creased along closed paths they are 
in a frustrated state [6 . Therefore, in order to balance 
their internal forces, they undergo buckling which leads 
to self-folding due to its highly constraint geometry. This 
suggests that mechanics should be fundamental to deter- 
mine the equilibrium configurations of these structures. 
Although Demaine et al [13] showed that a fold pattern of 
concentric squares cannot be rigidly folded without either 
adding additional folds or internally stretching the paper, 
the question of concentric circles is still open. Our numer- 
ical and analytic results strongly suggest that the conjec- 
ture of [13 that concentric circles can be folded without 
stretching is, indeed, true. A full proof, however, would 
need to show that the constraints of Eqs. (p!4|) and ([T5j) do 
not occur for sufficiently small but finite not just when 
w approaches zero. 
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